DFF 236/11/95 
November 10, 1995 



MAKING QUANTUM GRAVITY CALCULABLE 



R. Loll 1 

Sezione INFN di Firenze 
Largo E. Fermi 2 
1-50125 Firenze, Italy 

and 

Max- Planck- Institut fur Gravitationsphysik 
Schlaatzweg 1 
D- 14473 Potsdam, Germany 



Abstract 



We describe recent attempts at discretizing canonical quantum gravity in four dimensions 
in terms of a connection formulation. This includes a general introduction, a comparison 
between the real and complex connection approach, and a discussion of some open problems. 
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1. Introduction 



In this article I want to report on some recent developments in applying discretization 
methods to quantum gravity in four space-time dimensions. My emphasis will be on explain- 
ing some of the main ideas and motivations to the non-expert, rather than going much into 
the technical details, which for the most part can already be found in the literature. 

As is well known, there is as yet no complete consistent framework that would deserve 
to be called a theory of quantum gravity. However, there has been a great deal of activity 
during the last few years on the canonical quantization of gravity, and most of the material 
that follows will be related to these new developments. Having no very precise idea of what 
a complete theory of quantum gravity is eventually going to look like, it is not clear a priori 
what kind of quantities one is interested in calculating (or approximating) . One point of view 
is that the mathematical structure of this hypothetical theory would naturally lead to a set 
of preferred quantities that take on a particularly easy form, an idea borne out by recent 
results in the canonical approach to quantization. 

Our starting point will be the usual Einstein-Hilbert action for pure gravity, 



possibly with a cosmological constant term. In (1.1), is a non-degenerate, Lorentzian 
four-metric on the manifold M. The two basic approaches to the quantization of gravity, 
as given by the action (1.1), are the Lagrangian one using a path integral over all possible 
configurations of the system, and a Hamiltonian one, based on operator algebras of preferred 
functions or observables. In the first case, one is usually only able to treat four-metrics with 
Euclidean signature, which leaves one with the problem of having to continue the results to 
the Lorentzian regime. Although the two approaches may look rather different in concrete 
implementations, they nevertheless share some important fundamental questions in common. 

In the Lagrangian formulation, one works with functional integrals which, for example, 
are of the form 



in order to compute vacuum expectation values of functions O of the four- metric g M „, or 
rather of diffeomorphism equivalence classes [g^] of metrics, often called "geometries". The 
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integration in (1.2) is over some appropriate set of such geometries, with a suitable measure 
T>[g] that is supposed to depend only on equivalence classes of metrics, i.e., elements of the 
quotient space RiemM/DiffM. The usual strategy of quantum field theory would be to 
expand (1.2) perturbatively around the flat Minkowskian (or rather Euclidean) metric rj^, 
by considering metrics of the form g^ v = rj^ + h^ v that differ only by a small amount 
from the flat metric. However, for quantum gravity this strategy fails, since the theory turns 
out to be non-renormalizable. One therefore has to try to give meaning to the integral in 
some other, non-perturbative way. The crucial inputs in any such attempt are the space of 
quantum configurations to be integrated over, and the explicit form of the measure V [g]. 

In a Hamiltonian context, the central problem is to identify and quantize a preferred set 
of classical variables or observables. A typical such set consists of a sufficiently large Poisson 
algebra of functions which in the quantum theory is represented as an algebra of operators. 
As in the Lagrangian case, one has to determine the space of all quantum configurations 
and an inner product on this space in order to be able to compute scalar products and 
matrix elements of observables O. For example, in the metric approach, one seeks to make 
well-defined an expression like 



where for simplicity we have assumed that the four- manifold M was of the form M = ExS, 
with compact spatial slices S. The approaches to quantum gravity I want to discuss here 
take the viewpoint that one should start from discretized versions of quantities like (1.2) and 
(1.3), and then obtain the true theory in an appropriately taken (continuum) limit. 



2. Lagrangian approaches to discretization 

Two of the currently active research programs within the path integral approach are those 
of Regge calculus and dynamical triangulations. Quantum Regge calculus has developed 
from the ground-breaking idea of Regge [1] to do "general relativity without coordinates", 
by approximating a Riemannian manifold M by a piecewise linear, simplicial manifold A4. 
Since the manifold is piecewise flat, its curvature is concentrated on lower-dimensional hinges 
(in d = 4, on the two-dimensional triangles). A simplicial manifold Ai is characterized by an 
incidence matrix which describes the way in which individual simplices are glued together. 
In Regge calculus, one fixes an incidence matrix, i.e. a triangulation, and allows the edge- 
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lengths of the 1-simplices to vary. Edge-length configurations, subject to certain restrictions, 
then give rise to (discretized) Riemannian metrics. 



In the quantum theory, one tries to approximate the functional integral by the following 
product over all 1-simplices, 



where 5 Reggc is a discretized version of the Einstein-Hilbert action, depending on the link 
lengths U (see, for example, [2]). One central issue also here is the question of the right 
measure dfi(l). Since many edge length configurations correspond to the same metric, one 
expects the appearance of a non-trivial Jacobian J in dfi(l) = Y\ i dliJ{l). In two dimensions, 
progress has been reported recently in determining the measure factor analytically [3], how- 
ever, no similar results are known in higher dimensions. In numerical applications, the trivial 
measure Yli dk, or its scale-invariant version, Y\ i have been used, in the hope that they 
would lead to reasonable "effective measures" for the functional integral. 

The second approach, that of "dynamical triangulations" , may be regarded as a rigid 
version of Regge calculus. One starts from a set of equilateral simplices, with all edge lengths 
set to 1, say, and then considers all ways in which they may be glued together to obtain a 
manifold of a fixed topology (this implies restrictions on the possible glueings). The main 
observation then is that distinct such triangulations T give rise to different metric structures. 
One therefore introduces the functional integral 



where the Einstein-Hilbert Lagrangian (with cosmological constant) is discretized to S DT = 
c^N^T) — C2./V2(T). N2 and N4 denote the number of 2- and 4-simplices in the triangulation 
T, and the c» are bare coupling constants. This method was first applied in two dimensions, 
where a great deal of analytical results are available, and then generalized to higher dimen- 
sions. In this case, hardly any analytical results exist, and also the numerical analysis is 
much harder. The weight factor p(T) is usually set to 1. Relevant questions in the dynamical 
triangulations program are over which (inequivalent) triangulations the sum (2.2) should be 
performed, and whether in this way one obtains a uniform sample of Riemannian geometries 
(see [4] for recent reviews). 
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In either approach there are as yet no definite results about the existence of phase transi- 
tions and continuum limits in four dimensions, and therefore also the question of the restora- 
tion of the diffeomorphism symmetry (which is not present in the discretized formulation) in 
such a limit remains open. 



3. Discretizing in the canonical picture 

Let us now turn to the Hamiltonian approach, which will be the subject of the rest of 
this article. It will be indispensible for us to use a connection form A\ as our basic variable, 
instead of the usual three-metric g a b of the ADM approach. The basic pair of canonically 
conjugate variables will consist of an S'0(3)-connection A\{x) and a dreibein Ef . They are the 
exact analogues of the basic variables on the phase space of a Yang-Mills theory with gauge 
group G, where we would call the A^s gauge potentials and the Efs generalized electric 
fields. The far-reaching insight that also canonical gravity might be profitably described 
by gauge theory-like variables is due to Ashtekar [5]. The relation with the three- metric of 
the ADM canonical variable pair (g a b,TT ab ) is given by g ab = E^E bl . Note that, unlike in 
Yang-Mills theory, the Einstein-Hilbert action cannot be rewritten functional S[^A] of 
a four-dimensional connection only. The gauge-theoretic analogy therefore holds only after 
the 3+1-decomposition. 

In this formulation, and at a kinematical level, canonical gravity looks like a Yang-Mills 
theory with additional constraints, coming from the diffeomorphism symmetry not present 
in a gauge theory. Although the basic variables carry an internal index i, associated with 
the local S'0(3)-symmetry, physical quantities, like in Yang-Mills theory, have to be gauge- 
invariant. Denoting by A the space of all connections on E, this means that all physically 
relevant information is already contained in the quotient space A/G, the space of gauge 
equivalence classes of connections. 

Another input from Yang-Mills theory is the use of functionals of A labelled by closed 
curves 7 in the manifold S. These are the so-called Wilson loops, the traces of path-ordered 
exponentials of a connection A along a curve 7. Their importance comes from the fact 
that they are explicitly gauge-invariant and form an (over-complete) set of variables on the 
quotient space A/G- Their non-local nature, that may be regarded as a defect from a gauge- 
theoretical point of view, turns out to be an asset in gravity. This happens because spatial 
loops are extended objects with diffeomorphism-invariant properties: they can link, knot and 
intersect. 
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For our purposes, Wilson loops are interesting objects because they are readily discretized 
in a lattice framework. Consider a cubic, three-dimensional N x N x iV-lattice (with periodic 
boundary conditions), made up of vertices and one-dimensional, oriented links or edges l{. The 
link holonomy Ui of a link I can be thought of as the exponentiated integral of the connection 
A along that link. Since the A's are non-commuting matrices in the internal space, one has to 
take the path-ordered exponential of A. Note that no metric structure is necessary to perform 
the integration since A is a spatial one-form. We write Ui(A) = Pexp J l A, where one has 
to keep in mind that Ui takes values in the gauge group 5*0(3) and still has a non-trivial 
transformation behaviour under gauge transformations at the beginning and end point of the 
link I. To get rid of this dependence, one may combine a set of oriented lattice links so that 
they form a closed loop 7 on the lattice. Taking the trace of the resulting object, one obtains 
the lattice Wilson loop Tr U 1 := Tr \J\JJi 2 ■ ■ ■ Ui k . 



4. Real vs. complex connection formulation 

We do not want to derive the Ashtekar variables from scratch, but rather introduce them 
in a form that exhibits some special features. Starting from the Arnowitt-Deser-Misner vari- 
ables (g a i,,Tt ab ), one can derive an equivalent version of canonical gravity based on a pair 
(P?,Kl), with an internal S'0(3)-index i, where P" is again a dreibein variable, and the 
canonically conjugate momentum K l a is closely related to the intrinsic curvature of S. Com- 
pared with the usual ADM formulation, one obtains three additional first-class constraints, 
corresponding to the local 50(3)-symmetry. 

Consider now the canonical transformation 



where V is the S'0(3)-connection compatible with the triad P". The Ashtekar variables A l a 
and Ef obey the classical commutation relations 



since T l a is a function of Ef only. In this new canonical formulation, the first-class constraints 
read (up to constant factors) 




1 =. 
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Gi := V a E? = 0, (4.3a) 
H a := Fl h E\ = 0, (4.36) 
ft = -C,{3\^ k ~EfE)F ahk + 2(/3 2 C - l)E^(Ai - r»)(^ - rj), (4.3c) 



where Gi are the three Gauss law constraints (V = V(A) is the covariant derivative with 
respect to A), the H a are the three spatial diffeomorphism constraints (F l ab = d a A\ — d\,A % a + 
e 1 jkA J a Afr is the field strength of the connection A), and H is the Hamiltonian constraint. 
The algebraic expression for H contains two free parameters, (3 and There are two possible 
values for (, where 



_ j +1 44> Euclidean signature; 
1—1 44> Lorentzian signature. 

For special choices of f3 one can drastically simplify the form of the Hamiltonian H. In 
fact, if (/3 2 £ — 1) = 0, the second term in (4.3c) drops out and one is left with a polynomial 
Hamiltonian function. This leads to the following choices for (3: 



(R 2 t i) — o <^> / @ = = ' = ^' Euclidean signature; 
^ ' I f3 = ±i, Lorentzian signature. 

Our primary interest lies of course with the Lorentzian case. Nevertheless it is interesting 
to observe how the difference between a Lorentzian and a Euclidean signature for the space- 
time metric is reflected in the canonical formalism. The derivation for H in this compact 
form is a slightly modified version of the one given by Barbero [6]. It follows that if (in the 
Lorentzian case) we want to simplify the Hamiltonian constraint, we have to use a complex 
canonical transformation (i.e. (4.1) with (3 = ±i). This is the choice originally made by 
Ashtekar [5], and it leads to a complex »SO(3)-connection A\ as one of the basic variables. 
This gives rise to complications in the sense that one has to keep track of a set of reality 
conditions, to ensure that real, and not complex gravity is described. On the bonus side, 
all of the constraints are now polynomial in the basic variables A and E. This is important 
in the quantization, because polynomiality of a classical phase space function simplifies the 
search for a corresponding quantum operator. (Miraculously, if the real world were described 
by Euclidean metrics, one could achieve this simplification without ever leaving the real 
domain.) 
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The alternative that emerges from this picture (for £ = —1) is to set (3 = —1, say, so 
that A l a (x) is a real SO (3)-connection, and the basic Poisson bracket relations (4.2) have the 
standard form. As shown in [6], the Hamiltonian can be rewritten as 

~H = e^EfE* (F abk (A) - 2R abk (T)) = 0, (4.4) 

where R is the curvature of the connection T, R a b l = d a T % b — c^r* + e^r^r^. Because of the 
T-dependence, H is non-polynomial in the basic variables. It is rather similar in structure to 
the usual ADM Hamiltonian, 

H ADM [g,n] = -^g^R + ±^~ ab _ 1~ 2) = 0> (45) 

in that it can be brought into an equivalent form which is polynomial modulo powers of the 
square root of the determinant of the metric, ^fg = \J \E\. However, as is well known, its 
non-polynomiality is one important reason that makes the ADM Hamiltonian _£/~ ADM so hard 
to quantize, and has lead to an impasse in the canonical quantization in the metric approach. 
So it seems as if nothing had been gained by using this real connection formulation. We will 
be returning to this issue in due course. 

Let us summarize the two basic alternatives for the connection formulation of Lorentzian 
gravity: either one works with a simple, polynomial Hamiltonian, leading to complex con- 
nection variables and the need for reality conditions, or one stays within the real framework, 
but has to live with a more complicated, non-polynomial Hamiltonian constraint. 



5. The lattice discretization 

Having set the stage for the connection formulation a la Ashtekar in the continuum, we 
now come back to the discretized lattice formulation. The basic ingredients are those of the 
customary Kogut-Susskind approach of Hamiltonian lattice gauge theory [7]. For simplicity, 
we will work with the double cover ST/ (2) of SO(3), which in any case is necessary when 
considering matter coupling. For the complex case, the gauge group should of course be 
SC/(2) C = 51/(2, C). The lattice analogues of the variables (A % a ,Ef) are the link holonomies 
Ui introduced earlier, and the link momenta Pi(l). We label vertices of the cubic lattice by a 
triplet n = (ni, n^-, n^), = 0, 1, . . . , iV = 0, and links / by their origin n and a direction a. 
With this notation, the algebra of the basic quantum operators is given by 
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[U A B (n,a),Uc D (m,b)} =0 

[p i (n,a),U A c (m,b)} = -^5 nm 5 di) T iA B U B c (n,a) (5.1) 
[pi(n,a),p j (m,b)} = i S nm S ai e ijk p k (n, a), 

where the indices A,B,... are those of the defining representation in terms of 2 x 2-matrices. 
The commutators (5.1) are the exact quantization of the corresponding classical Poisson 
bracket relations. The holonomy operators Ua b act by multiplication, UA B ip(g) = UA B ip(g), 
and the link momenta are represented by the differential operators pi = —iTiA B UB C (^j)c A ■ 
The Tj are the usual Pauli c-matrices, rescaled by a factor i, and g stands for a group element 
(not to be confused with the same symbol for the metric used elsewhere in this article). 

This is the right moment to discuss some of the special features of lattice gravity, as 
opposed to lattice gauge theory. Firstly, in the real connection formulation, the total lattice 
Hilbert space H is given by the direct product over all links of the square-integrable func- 
tions on the group, H = ®iL 2 (SU(2), dg), with dg the Haar measure on SU(2). However, in 
the complex 5L(2, (D)-case, we immediately encounter the problem that there is no analogue 
of the left- and right-invariant Haar measure, because the group is non-compact. Secondly, 
we have already mentioned that gravity possesses an extra symmetry, the spatial diffeomor- 
phisms. This symmetry is badly "broken" by the lattice discretization (as it is by any other 
discretization), i.e. it is no longer present. Usually this leads one to the requirement that the 
symmetry should be restored in the continuum limit, when all link lengths shrink uniformly 
to zero, and/or the number of vertices goes to infinity. 

Until recently, no possible resolutions to the scalar product problem in the complex case 
were known. (The real case has so far not been considered.) Some investigations of lattice 
versions of Ashtekar gravity were made in [8-10]. Generally speaking, the lattice computations 
are rather lengthy, and not many results have been obtained concerning the algebra of lattice 
constraints, both classically and quantum-mechanically (in the latter case one also has to 
check for possible anomalies), and the continuum limit. Renteln was able to show that a 
particular discretized version of the quantum algebra of spatial diffeomorphism constraints 
reproduces the continuum algebra [9]. 

What comes to our aid in defining a scalar product for the SX(2, (D)-case is a construction 
by Hall [11], namely, an integral transform 
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C t :L 2 (SU{2),dg) - L 2 (SL(2, <D), du t ) ho1 

(5.2) 



[C t (/)](5c) == / dgf(g)p t (g-\ 
JSU(2) 



'St/(2) 

where the target space consists only of the holomorphic square-integrable functions on the 
complexified group ((5.2) is to be thought of as a non-linear analogue of the Segal-Bargmann 
transform L 2 (JR n ,dx) — > L 2 (C n , d/i) ho1 , where c#// is the Gaussian measure). The measure 
in the holomorphic representation is the heat kernel measure, depending on a positive, real 
parameter t, and p t is the heat kernel itself. 

These results have been used in [12] to set up a holomorphic lattice representation. 
Since the Hall transform is an isomorphism of Hilbert spaces, it induces a scalar product on 
SL(2, (D)-lattice wave functions (an explicit coordinate expression for dv t is not known), and 
it turns out that Wilson loop states are indeed square-integrable with respect to this inner 
product. The next step then is to find an appropriate lattice discretization H dlscr and look 
for solutions to the Wheeler-DeWitt equation on the lattice, H^^ip^j) = 0, and observables 
O satisfying [.ff" dlscr , O] = 0. By ^(7) we denote a generic linear combination of Wilson loop 
states on the lattice. If one thinks of taking eventually a continuum limit by letting the 
lattice spacing a go to zero, one has to solve these equations only to lowest order in a. In 
spite of the finite dimensionality of the lattice discretization, it is rather non-trivial to make 
progress on these questions. 

What is still needed in this context is a systematic investigation of all possible discretized 
Hamiltonians and scalar products. The choice of the classical H dlscr is non-unique because 
the criterion that H dlsCT ^5 H cont in the limit does not fix H dlsCT uniquely. Additional 
ambiguities arise through different operator orderings of the quantum Hamiltonian H dlscr . 
As for the scalar product, there are indications that it is too simple to reflect in a direct way 
the reality conditions A + = 2Y(E) of the classical phase space formulation. Rather, it 
corresponds to reality conditions of the form = A + 2E. The implicit metric dependence 
of this latter expression could lead to problems with the continuum limit in lattice gravity. 
Such a dependence is also present in the holomorphic continuum formulation [13], where one 
has to introduce an auxiliary function that depends on the lengths (!) of loop segments. Still 
this does not necessarily mean that the holomorphic Hilbert space cannot be put to some use 
in complex Ashtekar gravity. 

For specific choices of the discretized Hamiltonian, some solutions to the Wheeler-DeWitt 
equation have been found. The simplest ones depend on so-called Polyakov loops [12], which 
are straight closed lattice loops winding around the lattice in one of the three directions (this 
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is possible because of periodicity). More recently, other solutions have been constructed [14], 
which are labelled by multiple plaquette loops (a plaquette on the lattice is a smallest loop, 
made up of four links). Although they solve the Wheeler-DeWitt equation, they are not 
square-integrable with respect to the holomorphic scalar product. There may of course exist 
other scalar products in which those solutions have finite norm. 



6. Computing volumes 

Apart from the issue of the scalar product, there is another reason why the solutions 
to the discretized Wheeler-DeWitt equation found so far are probably not very interesting 
physically. In order to understand this, we need to study the so-called volume operator. It 
comes from the classical volume function, 



V(TZ) = J^d 3 x = J^d 3 x ^\e abc & k E?Effi\, (6.1) 

which measures the volume of a spatial region 1Z C S. Note that V is not an observable in 
pure gravity in the sense that it does not Poisson-commute with the Hamiltonian constraint 
H) . It however becomes an observable if the region 1Z is defined intrinsically in some way by 
the presence of matter fields [15]. This quantity can be studied also on the lattice, because 
it can be discretized straightforwardly and written as a function of the lattice momenta, 



v iatt = J^\e abc e^ k Pi (n,a) Pj (n,b) Pk (n,c)\, (6.2) 

where the integral has been substituted by a sum over all vertices lying within a lattice region 
7Z. Because of the presence of both the modulus and the square root, it is not a priori obvious 
how to quantize (6.2). However, in this case we are lucky, because both in the real and the 
holomorphic representation the momenta p are selfadjoint operators [16]. One can therefore 
go to a basis where the terms cubic in p under the modulus are diagonal operators, and the 
modulus and square root can be defined in terms of their eigenvalues. 

In order to diagonalize the volume operator, it is very convenient to use a particular set of 
gauge-invariant states in the loop representation, the (generalized) spin network states, as was 
first pointed out in [15]. A spin network state may be represented as a certain antisymmetrized 
linear combination of Wilson loop states. In order to describe a spin network on the lattice, 
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one needs to associate an "occupation number" G 2 + with each lattice link (subject to 
a number of restrictions). This can be represented by drawing lines through the link I. 
To obtain a multiple Wilson loop, one has to connect these lines at the lattice vertices in 
a gauge-invariant manner. The contractors at the vertices are the second ingredient in the 
definition of a spin network state. Details of this construction can be found in [15-17]. 

It can be shown that the spin network states span the space of all gauge-invariant 
lattice functions [18,17]. The operator V latt is "almost diagonal" in terms of the spin 
networks, i.e. the diagonalization is reduced to diagonalizations within finite-dimensional 
eigenspaces of V latt ; and its spectrum is discrete [15]. In determining the eigenvalues of 
the volume operator explicitly, it suffices to study the action of the operator D(n) := 
e abc e l i k pi(n, a)pj(n, b)pk(n, c) around a vertex n. Moreover it is useful to consider separately 
cases of different valence, where the valence counts the number of links with non-vanishing 

meeting at n. 

Because of the antisymmetry of the operator D(n), it annihilates spin networks which 
are two-valent at n. Although it is not immediately obvious, one can construct a general 
argument that also trivalent vertices are annihilated, i.e. spin network states that possess at 
most trivalent intersections have "zero volume" [16]. This is the reason why we considered 
the solutions for the Hamiltonian constraint equation described above as probably uninterest- 
ing. One obtains the first non-trivial eigenvalues for four-valent intersections. Some explicit 
spectral formulae for such spin networks were derived in [19]. One technical problem one 
encounters in the diagonalization is the occurrence of linearly dependent spin network states 
(because not all ways of contracting flux lines at vertices are independent). These states have 
to be eliminated to avoid spurious eigenvalues. 



7. Outlook, and a new perspective 

Apart from yielding interesting information about the volume operator itself, the con- 
struction described in the previous section also throws new light on the real connection 
formulation. We mentioned earlier that Barbero's Hamiltonian can be written in polynomial 
form up to powers of y/g. Now observe that in the diagonalized spin network basis, where 
the volume operator is diagonal, also the determinant g (which may be written as a cubic 
polynomial in the momenta) is diagonal. This means that in this basis it is straightforward 
to define all operators of the kind x /g k . It turns out that this enables one to quantize a 
discretized version of the Hamiltonian (4.4) on the lattice [20]. 
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Of course this quantum Hamiltonian is still more complicated than the one of the complex 
Ashtekar formulation, H = e l: > k E'j L EjF a bk, in that it has an additional potential term. On the 
other hand, there are no problems with the choice of a scalar product which is simply given 
by the ordinary Haar measure on SU(2), as in lattice Yang-Mills theory. However, contrary 
to the ADM formulation, there is now an obvious way to quantize arbitrary functions of the 
determinant of the metric. 

If one wants to pursue this line of thought further, and try to solve the Wheeler-DeWitt 
equation, say, it is clear that one a) needs an efficient general method for generating volume 
eigenstates, and b) must eliminate the zero-volume states because of the inverse powers of 
the metric occurring in the Barbero Hamiltonian. Thus, unlike in complex Ashtekar gravity, 
in the real context one cannot extrapolate smoothly to the case where the (quantized) three- 
metric is allowed to be degenerate. 

To summarize, there has been progress in Hamiltonian lattice gravity, although many 
details remain to be worked out. Crucial in all recent developments has been the use of 
connection variables, and the associated gauge-invariant Wilson loops. There are two basic 
alternatives to proceed in the lattice framework. The first one is that corresponding to 
complex A; there it remains to find solutions "with volume" to Hip = (although it is 
not strictly necessary to work in terms of spin network states), and to investigate different 
Hamiltonians and scalar products. The second way is that based on a real connection A] as we 
have seen, in this case the use of a diagonalized spin network basis is mandatory; the search 
for solutions to the Wheeler-DeWitt equation has only just begun [20]. Once some solutions 
have been found, it would be very interesting to compare the results of both the real and the 
complex formulation, and to see whether they are equivalent. Other issues that remain to 
be addressed are whether and how the continuum limit should be taken, and whether and 
how diffeomorphism invariance should be incorporated in the lattice picture. Also it would 
be most interesting to try to study matter couplings in the same framework. 

Acknowledgement. I am grateful to the organizers of the Sintra meeting for their hospitality, 
as well as to the members of the Max-Planck-Institut where these notes were written up. 
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